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Abstract: A new, modular toolkit for creating simulations of 2D X-ray pixel detectors, X-CSIT 
(X-ray Camera Simulation Toolkit), is being developed. The toolkit uses three sequential simula¬ 
tions of detector processes which model photon interactions, electron charge cloud spreading with 
a high charge density plasma model and common electronic components used in detector readout. 
In addition, because of the wide variety in pixel detector design, X-CSIT has been designed as a 
modular platform so that existing functions can be modified or additional functionality added if the 
specific design of a detector demands it. X-CSIT will be used to create simulations of the detectors 
at the European XFEL, including three bespoke 2D detectors: the Adaptive Gain Integrating Pixel 
Detector (AGIPD), Large Pixel Detector (LPD) and DePFET Sensor with Signal Compression 
(DSSC). These simulations will be used by the detector group at the European XFEL for detector 
characterisation and calibration. For this purpose, X-CSIT has been integrated into the European 
XFEL’s software framework, Karabo. This will further make it available to users to aid with the 
planning of experiments and analysis of data. In addition, X-CSIT will be released as a standalone, 
open source version for other users, collaborations and groups intending to create simulations of 
their own detectors. 

Keywords: XFEL; Simulation; PIXEL; Pixel detectors; Semiconductor detectors. X-ray 
detectors. 
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1. Introduction 

X-CSIT is a toolkit for creating simulations of 2D semiconductor pixel detectors. It provides 
physics simulations from incident photons to the readout electronics. The toolkit is modular and 
adaptable, as its design goal is to provide a common simulation framework for the wide variety 
of detectors to be used at the European XEEL (XEEL.EU) [U, 0, 0, 0], which are summarised in 
table |T|. Detectors to be simulated by X-CSIT at the European XEEL include: the Adaptive Gain 
Integrating Pixel Detector (AGIPD) 0, 0], the Large Pixel Detector (LPD) EL 0], the DePEET 
sensor with Signal Compression (DSSC) [0, [13, HID as well as pnCCDs [jI3, [HD and EastCCDs W\M . 

X-CSIT will be integrated into the European XEEL’s software framework Karabo flTTD . As 
part of the integration, XEEL.EU users will be provided with pre-built simulations of XEEL.EU 
detectors. Karabo will allow X-CSIT simulations to run on XEEL.EU’s computer network and 
transparently handle concurrent processing. Eurthermore, simulations can be tied to the data pro¬ 
cessing pipelines implemented in Karabo, allowing simulated data to be processed in the same way 
as measured data. 

2. Design of X-CSIT 

2.1 Design goals and objectives 

X-CSIT was conceived from the desire to have a single common, well validated simulation of the 
detectors to be used at European XEEL, see table ^ capable of handling large amounts of data in a 
streaming fashion. This requirement necessitates a versatile and adaptable simulation and contrasts 
with existing simulation tools, which usually require processing all data in sequential steps, e.g. 
HORUS [1T3I1 or MEGALib flTBl . X-CSIT achieves this by providing a set of physics simulations 
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Table 1: Selected specifications for some of the detectors to be used at the European XFEL. Note 
that the pnCCD and EastCCD can also be used for imaging spectroscopy. Accordingly, an energy 
resolution dE is given for these detectors. 



LPD 

AGIPD 

DSSC 

pnCCD 

FastCCD 

Pixels 

500/im 

square 

200/im 

square 

204;itm 

hexagonal 

75/im square 

30/im square 

Depth 

500/im 

500/im 

450;itm 

300/im 

200/im 

Dynamic 

range 

1x10^ at 12 
keV 

1x10^ at 12 
keV 

6000 at 1 keV 

1000 at 2 keV, 
dE: 130 eV at 

5.9 keV 

1000 at 500 eV, 
dE: 400 eV at 5 

keV 

Dynamic 
range tech¬ 
nique 

Triple gain 
profile 

Pre-amplifier 
chosen gain 

DePFET non¬ 
linear gain 

Linear 

Linear 

Sensor size 

32x128 pix¬ 
els 

512 X128 pix¬ 
els 

256 X128 pix¬ 
els 

200x128 pix¬ 
els 

1920x960 pix¬ 
els 

Photon 

energy range 

12 keV op¬ 
timal, 1-24 

keV 

3-13 keV 

0.5-6 keV op¬ 
timal, 0.5-24 
keV 

0.1-15 keV 

0.25-6 keV 


based upon user provided detector definitions. Because X-CSIT is based on an object orientated 
design, code can be customised without affecting the rest of the simulation. 

The inherent versatility of X-CSIT and its ability to create simulations of a wide variety of 
detectors makes it useful beyond its original scope of the detectors at XFEL.EU. Once testing has 
been performed, and the physics simulations used within X-CSIT have been validated, the code 
will be made available as open source. Other groups can benefit from X-CSIT, which will provide 
a faster and more convenient means of creating a simulation than starting from scratch, and will 
already be validated. 


2.2 Simulation details 

X-CSIT consists of three simulations (figure [T|), each covering a separate physical regime within 



Charge 

Simulation 


Plasma 

Search 


Charge 

Generation 




Charge 

Spreading 


Plasma 

Simulation 


/ 


Electronics 

Simulation 


/ // 


Noise 

\ Digitization 



Amplification^ 

1 Charge 
Transfer 


Figure 1: The layout and components of the X-CSIT sub-simulations. As input a list of primary 
particles to generate is passed to the particle simulation. After all simulation steps have been 
completed, images containing events in detector units are output. 
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the detector. Data passes sequentially through these simulations, although all three can be run 
concurrently on different data packets, thus constituting a pipeline. 

Between these simulations data is passed using container classes that define the interface that 
X-CSIT uses. By deriving from these classes, user-defined data storage classes may be defined. 
This allows for the expansion of the simulation, e.g. to include more detailed reporting of cer¬ 
tain simulation aspects, or storage of data in a user specified way. Karabo integration of X-CSIT 
at XFEL.EU utilises this functionality to make use of Karabo’s native data storage methods for 
optimisation of in-framework data transport. 

2.2.1 Particle simulation 

The particle simulation simulates individual energetic particles interacting with the experimental 
setup. To produce the presented results Geant4 [UBI, 13 version 10.0 was used for this simulation 
stage. Using the Livermore physics list fEOl . photon and energetic electron interactions including 
the photo electric effect, Compton and Rayleigh scattering as well as fluorescence and Auger emis¬ 
sions can be simulated. These processes have been validated for version 9.6 down to 1 keV EE E2]- 
In-house validation will be conducted for the Geant4 version integrated within X-CSIT. 

Users provide the particle simulation with a description of the physical geometry to be used 
by Geant4 and identify the sensitive detector regions. X-CSIT then manages the Geant4 simulation 
and handles the streaming of output data from Geant4 to the subsequent simulation components. 
In this way partial data can already be passed on to the charge and electronics simulations, rather 
than having to wait for the full Geant4 simulation to finish. 

2.2.2 Charge simulation 

The charge simulation models the behaviour of charge clouds in the semiconductor layer and the 
sharing of charge between pixels. Because of the large number of electron-hole pairs created, of 
the order of 275 per keV of initial energy in silicon, the simulation uses a statistical distribution, 
rather than a Monte Carlo simulation as used for the particle simulation. Additionally, two regimes 
are considered. 

In most cases separate charge clouds do not interact with each other noticeably, and even the 
influence of interactions between electrons and holes is small compared to that of drift and diffu¬ 
sion. If, due to the number of charges, only drift and diffusion need to be considered a statistical 
distribution of charge can be assumed. This distribution is a Gaussian normal distribution for which 
the standard deviation, < 7 , is proportional to the root of the distance between the interaction point 
and the collection point, d, over the electric field due to the bias voltage, E, as given by 



( 2 . 1 ) 


This is a valid approximation for strongly overbiased detectors, where interactions do not 
happen close to the electrode (d> 100 /rm) E3]- 

The two dimensional spread is calculated as a pair of orthogonal one dimensional Gaussian 
distributions in the X and Y axes. The proportion of charge crossing a pixel boundary can then be 
calculated with the Gaussian cumulative distribution (see figure 0). 
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Figure 2: Charge sharing calculations in X-CSIT. 


In the case of very high charge carrier densities the charges begin to screen themselves from 
the electric field inside the sensor, constituting a plasma. This plasma releases charge slower than 
in the non-plasma case, thus increasing the spread of the charge cloud as well as the charge col¬ 
lection time. Before the advent of FEL’s, charge clouds resulting from X-ray photon interactions 
in semiconductor detectors rarely reached the densities required for plasmas to occur lUM . and 
could frequently be neglected in simulations. In contrast, a simulation of plasma effects is required 
for use at XFEL.EU, since intensities reached there are expected to generate electron-hole plas¬ 
mas. The plasma model to be used will be developed empirically based on data due to be taken at 
XFEL.EU in 2015. 

2.2.3 Electronics simulation 

The electronics simulation consists of a set of modular components, each of which simulate an 
electronic component or effect such as amplification, digitisation or noise. A simulation of the 
electronics of a detector is created by chaining together these components to create a functional 
representation of the detector circuitry as seen by the collected charge moving through it. An 
example is shown in figure The empirical models used in these components require parameters 
given either as expected detector performance characteristics, or derived from calibration data taken 
from an existing detector. In the latter case either a specific detector or, by randomising the input 
data, a detector of a given detector type may be simulated. 

X-CSIT will come with a set of general components which can be arranged to model common 
detector front-end electronics designs. Additional components to model specific electronics can 
be added. At the European XFEL each component of the electronics simulation is integrated into 
Karabo as a separate device, conceptually similar to a plugin. In a future release of Karabo this will 
allow users to layout the electronics simulation of their detector graphically within Karabo’s GUI. 

3. Initial Testing 

Initial testing of X-CSIT has been conducted with data taken by pnSensor GmbH using a pnCCD \TA 
ITHl (for specifications see table [T|) and an Fe-55 source. A simulation of the pnCCD was imple- 
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(b) The proposed LPD electronics simulation compared to a diagram of the real 
electronics [0]. 


Figure 3: Electronics simulation diagrams for the pnCCD and LPD detectors. 
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Figure 4: Diagram of the setup of pnCCD Fe-55 measurements as simulated in X-CSIT. The alu¬ 
minum attenuator was used to reduce the photon intensity such that event pile-up on the sensor is 
avoided. 


mented in X-CSIT (see figure ^ and the data from the real and simulated detectors were run through 
the same analysis pipeline previously validated against an existing calibration and characterisation 
pipeline used for the CCD of the CERN Axion Solar Telescope fE5l1 . 

Figure ^ shows the histograms produced from the uncorrected measured and simulated data. 
The residuals between the two data sets shown are given by (measured — simulated) /\/measured. 
These spectra are not normalised with respect to each other. Instead, using the known source 
characteristics, the simulations were run to produce a number of primary photons that is equal 
to the number of photons incident on the detector in the measurements. Good agreement can be 
seen between the measured and simulated data above the noise peak (3000 ADU). The data in the 
regions just below the ^^Mn peaks are underestimated by the simulation, while the region below 
this, 8000 ADU to 9500 ADU, and the shoulder of the noise peak are overestimated. This could 
indicate that the simulation creates an excess of charge sharing events with the majority of the 
charge in the primary pixel of the split event. The low energy cutoff in the pnCCD is not simulated. 

The left panel of figure § shows a scatter plot of simulated events, where a single photon has 
shared charge between two pixels. The same features visible in the plot for the simulated data exist 
for the measured data (not shown), notably the Al fluorescence line, the escape peak, the ^^Mn 
Ka and Kp lines and pile-up events. The right panel shows the deviation between simulated and 
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Figure 5: Top panel: the uncorrected energy spectra for the measured and simulated data given 
in detector units. Lower panel: the relative deviation between the two, given in terms of 1 <7 
uncertainties. Visible peaks from left to right are: noise peak, A1 fluorescence, Kc^ and ^^Mn 
decay lines at 5.9keV and 6.5keV. 


Simulated 




Figure 6: Left panel: scatter plot giving the distribution of charge shared between the two pixels 
of a double pattern (simulated data). The data shown has been offset and common mode corrected. 
Points are colour coded by fractional Y position on the sensor, the position of the event along the 
readout direction. Right panel: Relative deviation between measured and simulated data given in 
terms of 1CJ uncertainties on the measured data. 


measured data in terms of 1 <7 uncertainties. In the central region the simulation slightly but system¬ 
atically underestimates the number of events by 2 < 7 , while at the edges a significant systematic 
overestimation occurs. This corresponds to the previous observation, that in the simulation a larger 
fraction of charge is attributed to the primary pixel of the split event. Additionally, it is apparent 
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that the simulation has fewer background events, i.e. events which cannot be attributed to the afore¬ 
mentioned features. Whilst only double patterns extending perpendicular to the readout direction 
are shown here, the behaviour is similar for patterns extending along the readout-direction. 

In figure 0 (a-c) spectra are shown for the different event types, characterised by the shape of 
the spread on the pixel grid. Additionally, first singles, i.e. singles that are also the first event in 
a readout column per transfer and thus has not been affected by the residual of a charge transfer 
of another event, are shown. These show qualitatively similar features on both the simulated and 
measured results in all four event types, however quantitatively the ratios of event types do not 
match. Events in three regions were compared for each event type: the low energy background 
region, ranging up to energies of 2000 ADU 1.58keV, the continuum region, in which split part¬ 
ners will largely be found (2000 ADU-6000 ADU ~ 4.75 keV) and the peak region, containing the 
^^Mn photo-peaks (> 6000 ADU). For each region mean deviations in terms of sigma uncertainties 
of the measured data have been calculated. These values are given in table 0. 

Overall, i.e. across the complete spectrum, triple and quadruple patterns are both over esti¬ 
mated (by 0.40% and 13.60%) while double patterns are underestimated (20.32%). The inclusion 
of background noise in the statistics of the measured singles and first singles but not the simu¬ 
lated singles makes comparison of these patterns (—51.03% and —52.73% in simulation) difficult. 
When considering only the immediate peak region for singles, i.e. events above 7000 ADU, thus 
eliminating the low energy background in the measured data, simulation singles and first singles 
are underestimated (—40.06% and —42.34%). 

For the individual regions it is found that the shape of the continuum is generally represented 
well by the simulation, although a consistent underestimation exists, except for quadruple split 
patterns. The peak region is consistently overestimated by the simulation, both for singles and 
higher multiplicity patterns. The large deviations in the background region can be attributed to 
the fact that a simplified geometry was simulated and as such scattered photons from experimental 
components other than the detector and the source are missing. This is a systematic effect, as is 
evident from the residuals shown in the plots and the large variance of deviations in this region. 

Figure 0 (d) shows a comparison between fully calibrated measured and simulated data, which 
has been corrected for offset, common-mode and charge transfer inefficiency (CTI), alongside fits 
to the ^^Mn and lines. The the numerical components of these fits and how they match to the 
expected emission lines are given in table 0. Additionally, the residuals between the two data sets 
are shown as (measured — simulated) /\/measured. The charge sharing excess in the simulation 
causes a deficit of singles, this causes higher uncertainties in the simulation fit. A match within Icr 
of fitted energy is seen between the simulated and measured data in both peaks as well as the width 
of the ^^Mn Kp line. The width of the fit of the ^^Mn line is larger for the simulated data set by 
a statistically significant amount, but the source of this deviation has not yet been identified. The 
deviation between the fully calibrated simulated and measured first singles match the deviations 
observed for non-calibrated first single events. 
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Figure 7: Plots (a-c): Energy histograms of the four different event patterns observed, defined 
by the shape of the spread on the pixel grid in any rotation, as shown in the legend at the top. 
Invalid corresponds to any other event shape, or patterns with a diagonally offset event adjacent to 
it. The data has been offset and common mode corrected. Each pixel of an event contributes to the 
histograms. Plot (d): Photo peak region showing a comparison of fully corrected data alongside 
the fit models used to determine the gain. 
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Table 2: Mean deviations between simulated and measured data in terms of 1 <7 uncertainties of the 
measured data for different event types and regions. The values in brackets give the variance of the 
data in a specific region. 


Region 

Background 

Continuum 

Peak 

Singles 

9.3 ( 80.1) 

1.3( 1.0) 

4.9(23.9) 

First singles 

14.1(191.4) 

1.6 ( 2.4) 

9.1(61.0) 

Doubles 

2.3 ( 33.2) 

3.5 ( 2.0) 

0.9 (14.4) 

Triples 

-2.1(190.5) 

2.7( 1.3) 

-4.2(28.3) 

Quads 

-14.8(424.2) 

-7.7 ( 8.6) 

-6.1(59.6) 

Invalids 

-11.1(958.8) 

-0.7 ( 2.4) 

-4.2(45.2) 


Table 3: Tabulated data of the fits to the ^^Mn peaks in figure ^ (d). 



Simulated 

Measured 

E (^^Mn Ka) 

(5892.78 ± 2.61) eV 

(5891.41 ± 2.46) eV 

FWHM (^^Mn K^) 

(151.63±0.58)eV 

(141.20±0.35)eV 

Ei^^Mn K^) 

(6490.19 ± 3.37) eV 

(6488.13 ± 3.05) eV 

FWHM (^^Mn K^) 

(158.50±1.81)eV 

(157.67±1.17)eV 


4. Conclusion and outlook 

X-CSIT is a toolkit for creating simulations of 2D semiconductor pixel detectors. These simula¬ 
tions include the simulation of photon interactions in the semi-conductor material, charge sharing 
between pixels and the response of the electronic readout. X-CSIT has been designed for adapt¬ 
ability to different detector designs. While this is required to take into account the variations in 
design of the European XFEL detectors it makes the toolkit useful for other groups as well. At the 
European XEEL the integration of X-CSIT into Karabo will make it available to XFEL.EU users 
with pre-set up simulations of the available detectors. 

Currently, an early working version of X-CSIT has been used to simulate a pnCCD which is 
also being used to do measurements for early testing and comparisons. Results indicate reasonable 
agreement between measurements and simulation, within 3(7 on the uncorrected spectrum. Sys¬ 
tematic deviations, especially concerning charge sharing, remain and are under investigation. Later, 
as more sources and detectors become available at European XFEL, additional physics elements of 
the toolkit will be validated. After the toolkit has been validated it will be made available as open 
source for other groups to use. 
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